Aharonov-Bohm oscillations in disordered topological insulator nanowires 
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A direct signature of electron transport at the metallic surface of a topological insulator is the 
Aharonov-Bohm oscillation observed in a recent study of Bi2Se3 nanowires [Peng et al., Nature 
Mater. 9, 225 (2010)] where conductance was found to oscillate as a function of magnetic flux 
through the wire, with a period of one flux quantum cftQ = h/e and maximum conductance at zero 
flux. This seemingly agrees neither with diffusive theory, which would predict a period of half a 
flux quantum, nor with ballistic theory, which in the simplest form predicts a period of 0o but 
a minimum at zero flux due to a nontrivial Berry phase in topological insulators. We show how 
h/e and h/2e flux oscillations of the conductance depend on doping and disorder strength, provide 
a possible explanation for the experiments, and discuss further experiments that could verify the 
theory. 

PACS numbers: 73.25. +i, 73.20.Fz, 73.43. Qt 



The characteristic feature of a strong 3D topological in- 
sulator (TI) is the presence of a conducting surface that 
is topologically protected from Anderson localization by 
time-reversal-invariant disorder [1]. In general the sur- 
face state has an odd number of Dirac points in the en- 
ergy spectrum, with the simplest case of a single Dirac 
point being realized at the (111) surface of Bi2Se3 [51[3]. 
While the presence of this surface metallic state has been 
demonstrated convincingly using surface probes, notably 
angle-resolved photoemission spectroscopy [2|, studies of 
the transport properties of these surfaces are rare HHH] ■ 
Topological insulator nanowire with perfectly insulating 
bulk realizes an ideal hollow metallic cylinder with a di- 
ameter large enough that it is easy to thread a large mag- 
netic flux through its core [6] . The magnetoconductance 
of such wires not only reflects the fundamental effects of 
normal metal physics such as the Aharonov-Bohm effect 
and weak localization, but can also indirectly probe the 
existence of a nontrivial Berry phase. 

The magnetic flux affects the transport properties of 
the metal surface through the Aharonov-Bohm effect: the 
wave function of the particles picks up a phase of 27T(p/(j)o 
going around the circumference, with (f> the total flux 
through the cylinder and 00 = h/e the flux quantum [5]. 
There are two inequivalent values of flux that do not 
break time-reversal symmetry in the surface: and 0o/2 
(up to integer multiples of c/jq). In normal metals there 
is no fundamental difference between the two values, but 
in TI nanowires there is: only one of these values allows 
for a state at the Dirac point (cf. Figjl]). When there is 
a state at the Dirac point, the total number of modes is 
odd and time-reversal symmetry requires the presence of 
a single perfectly transmitted mode [10] contributing con- 
ductance e^ /h. If the contribution of all other modes to 





Figure 1: A schematic of the band structure of a topologi- 
cal insulator nanowire in the presence of a magnetic flux (j). 
Because of a nontrivial Berry phase of the particle a gapless 
mode is only realized at a flux with half integral flux quanta. 

the conductance is exponentially suppressed, which hap- 
pens for example at the Dirac point for ballistic wires 
with a small aspect ratio (length ^ circumference), the 
conductance is dominated by the presence or absence of 
the perfectly transmitted mode. In a fiat surface with a 
Dirac point, such as graphene, the zero mode is realized 
at zero flux, while in a curved surface as in TI nanowires, 
it is realized at (j)o/2 pTMl3) . This phase shift occurs 
because the particle spin is constrained to lie in the tan- 
gent plane to the surface and thus a particle picks up a 
Berry phase of tt due to the 27r rotation of the spin as it 
goes around the surface [H]. The magnetoconductance 
of an undoped ballistic TI nanowire is thus expected to 
oscillate with a period of (j)Q and a maximum at 4> = </'o/2. 
The situation is different in the presence of disorder 
strong enough that transport is diffusive. The classical 
conductance acquires a quantum correction due to in- 
terference between time reversed paths. The phase dif- 
ference between a path going once clockwise around the 
cylinder and one going anticlockwise is 2 x 2t:(J)/(J3q and 



the conductance thus oscillates with a period of (?io/2 [5]. 
Whether the conductance has a maximum or a minimum 
at zero flux depends on the presence (weak antilocaliza- 
tion) or absence (weak localization) of spin-orbit cou- 
pling. By their very nature, topological insulators have 
strong spin-orbit coupling and thus show weak antilocal- 
ization with a maximum conductance at zero flux. 

Surprisingly, however, in the recent transport experi- 
ments on Bi2Se3 nanowires neither of the above scenarios 
seems to be realized [5]. In the experiment the weak an- 
tilocalization induced 4>q/2 period is essentially absent, 
while the 0o periodicity is clearly seen. However, the 
conductance has a maximum at = rather than at the 
expected <j) = 4'ol'^- Although none of the available TI 
are particularly good bulk insulators, an explanation of 
this discrepancy from bulk properties is unlikely, since 
there is no compelling reason for the bulk conductance 
to show sharp flux periodicity determined by the cross- 
sectional area. 

In this paper we provide a theoretical study of the 
transport properties of the surface in the presence of par- 
allel flux and time-reversal preserving disorder, combin- 
ing analytical estimates with numerical simulation. By 
studying a pure surface theory our results are not com- 
plicated by bulk contributions, making it easier to dis- 
entangle the surface and bulk properties in the exper- 
iments. We show that in the regime of weak disorder 
and nonzero doping our theory is consistent with the ex- 
perimental result. We begin by briefly explaining our 
theoretical model and calculations before presenting our 
results and discussing them in relation to experiments. 

The transport properties of the surface are determined 
by the Dirac equation 



[wp • a + V[y)\ i/j = eip. 



(1) 



"■ — (o'a;,crj,) are the Pauli sigma matrices and v is the 
Fermi velocity. The Fermi energy e is determined by the 
density of surface charge carriers and can in principle be 
tuned by doping. Vir) is the disorder potential which we 
take to be Gaussian correlated 
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(2) 



with correlation length ^. The exact form of the corre- 
lator is not important in obtaining our results. Kq is 
the dimensionless measure of the disorder strength. We 
take < a; < L to be the coordinate along the wire 
and < y < W^ as the circumferential coordinate. The 
magnetic flux </> is absorbed into the boundary condition 



^p{x, y + W)^ V(a;, y)e*(2'^^/'*"^ 



(3) 



The extra factor of tt is the curvature-induced Berry 
phase. Despite the large g factor the Zeeman coupling 
of the magnetic fleld and the spin is not expected to be 
of relevance. Since the field is parallel to the surface this 
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Figure 2: Conductance versus Fermi energy for three different 
disorder strengths Kt^ and for three representative values of 
magnetic flux. Here W = 100^ and L — 200^. At low disorder 
(K=0.2), whether the conductance at = or = (j}o/2 is 
larger is highly sensitive to the location of the Fermi level. 
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Figure 3: Conductance at the Dirac point e = as a function 
of magnetic flux with W = 100$, L/\ = 100 (sohd hues), 150 
(dashed hnes) and 200 (dotted lines), and two values for the 
disorder strength. 

simply shifts the band structure in Fig. [l] unlike a nor- 
mal Zeeman fleld, which opens up a gap. We therefore 
ignore Zeeman coupling in this work 

The scattering matrix is obtained using the transfer 
matrix method of Ref. 1151 which in turn gives the con- 
ductance G through the Landauer formula. The total 
transfer matrix T relates the wave function on the left 
{x = 0) to the wave function on the right {x = L) and 
is obtained as a product of N transfer matrices 5Tj that 
propagate the wave function from x — {j ~ 1)L/N to 
X = jL/N with j — 1, . . . ,N. The STj are obtained by 
solving the Dirac equation (fTl) . The dimension of 6Tj is 
determined by the number 2M -I- 1 of transverse modes 
(j„ = 27r(n + (f)/(j)o + 1/2) /W, n = -M, M -h 1, ..., Af in- 
cluded in the calculation. We take A^ and M large enough 
that the conductance no longer depends on them. For ex- 
plicit expressions for the matrices STj and further details 
we refer to Ref. 15, We average over ~ 10"^ disorder re- 



alizations. 

The model (fTl) has been studied extensively in the limit 
of large aspect ratio W/L >• 1, in the context of graphene 
in the presence of valley-preserving scalar disorder (for a 
recent review see Ref . \W\\ . In this limit the conductance 
is independent of boundary condition, and therefore of 
flux, and is topologically protected from Anderson local- 
ization, growing logarithmically with system size indica- 
tive of weak antilocalization [T51 [T7] . The opposite limit 
of very small aspect ratio W/L ^ 1 describes carbon 
nanotubes if in addition W is very small (smaller than 
the mean free path). Carbon nanotubes are fundamen- 
tally different than TI wires in that the Berry phase is 
absent. Also, due to their small diameter it is challenging 
to thread through them a flux of the order of or larger 
than a flux quantum. For a review of related theoreti- 
cal and experimental work see Ref. [T51 Here we study 
the intermediate regime of aspect ratio W/L ^ 1, which 
is the regime relevant to the experiment of Ref. 6, In 
this regime the conductance does depend on the flux and 
this dependence in the presence of disorder has to our 
knowledge not been studied systematically before. 

We now present and interpret the results of our numer- 
ical simulation. In Fig. [2] we plot the density dependence 
of the conductance for a few different disorder strengths 
and a few values for the magnetic flux. Three qualita- 
tively different regimes are observed: (i) The small dop- 
ing regime close to the Dirac point, and large doping 
regime with (ii) weak disorder (small Kq) and (iii) large 
disorder (large Kq). 

In the small doping regime, the conductance is gener- 
ally small compared to the conductance quantum e^//i, 
and the physics in this regime is thus dominated by the 
possibility of a perfectly transmitted mode. In Fig. Islwe 
show the flux dependence of the conductance for a fixed 
diameter and a few wire lengths. Because of the Berry 
phase the conductance peaks at (jj = 0o/2 with conduc- 
tance of about e^/ft.. For other values of flux the conduc- 
tance is exponentially suppressed (with aspect ratio) and 
goes to zero in the limit of W/L — > 0. Interestingly, the 
conductance increases with increasing disorder strength 
for a fixed system size, at the same time going to zero 
with increasing system size. 

For weak disorder, away from the Dirac point, the mag- 
netoconductance generically oscillates with a period of 
(f)Q with an amplitude that depends on e. The maximum 
conductance is either at (j) = or at cf) — (t>Q/2 depending 
on the amount of doping. (Without doping, the maxi- 
mum is at = 00 /2, as discussed in the introduction.) 
In the inset to Fig. Ill we plot the amplitudes (5Gi/2 and 
SGi of the 00/2 ana0o-periodic oscillations of the mag- 
netoconductance, 

.^ G(0 ^ 0o/2) + G(0 = 0) - 2G(0 = 0o/4) 

SGi = , (4a) 




Figure 4: Main panel: Comparison of the analytical expres- 
sion ([6| for the weak antilocalization correction as a function 
of aspect ratio, with numerical results. The numerical data 
are obtained for a fixed W = 40^, Kq — 1 and e^/vh = 3. 
Inset: The conductance amplitudes SG\/2 and SGi of Eq. Q 
as a function of disorder strength Kq at e^/vh = 0.308 and 
for L — 2W = 200^. The structures at small Ko can be un- 
derstood as Fabry-Perot resonances. As disorder increases, 
the half-flux oscillations dominate the one-flux oscillations. 

as a function of disorder strength for e = 0.308/iw/^ as 
an example of a doping level where the conductance has 
a maximum at = (as in the experiment [6]). 5Ga 
is a measure of the strength of the a0o period in the 
magnetoconductance. 

For very small values of Kq the conductance can have 
sharp Fabry-Perot resonances that leads to the compli- 
cated nonmonotonic behavior seen in Fig. |4] The reso- 
nances are very sensitive to disorder and disappear al- 
ready at relatively small values of Kq [TU] . With increas- 
ing disorder strength the nonmonotonic dependence on 
the density is smoothed out and G(0 — 0o/2) becomes 
equal to G(0 — 0) and the half-flux quantum period 
starts to dominate over the one-flux quantum period. 
Away from these values of the flux, time-reversal sym- 
metry is broken and the weak antilocalization correction 
to the conductance gives a 0o/2 period of the magneto- 
conductance. 

We can estimate the disorder strength at which the 
crossover between the 0o and 0o/2 period happens with 
the following argument. The elastic linewidth broadening 
is given by [20] 
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<5Gi = G(0 = 0) - G(0 = 0/2), 



(4b) 



with kp = e/{hv). When this broadening is of the order 
of the mean level spacing 2hv7T/W, the oscillations in the 
density of states that cause the one-flux quantum period 
are washed out. For the parameters in Fig. |4]this gives 
a crossover value Kc ~ 0.2 in reasonable agreement with 
the data shown in Fig. [4J 



In the limit of short and wide wire the conductance is 
independent of the flux, and 5Gi/2 goes to zero. Intu- 
itively this is because the particles leave the wire before 
they have time to go around the circumference of the 
wire. Adapting the calculation of the weak antilocaliza- 
tion correction [|21. to the boundary conditions appropri- 
ate to the cylinder geometry one finds 



5G 






(6) 



with 



oa — — log — H — > cos 
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up to an L- and (/)-independent constant. In Fig. |4] we 
compare the analytical expression m^ to the numerical 
calculation with good agreement. 

With a good understanding of the different transport 
regimes we are now in a position to discuss the findings 
of the experiment of Peng et al. [6j in the context of our 
results. A magnetoconductance with a period of 4>o and a 
maximum at zero flux is realized in our model at doping 
large enough that the conductance is larger than e^ /h 
and in the presence of weak disorder such that Fabry- 
Perot resonances are washed out but the dynamics are 
not yet fully diffusive. It is reasonable to expect the 
samples to be doped as there is no reason for the Fermi 
energy to coincide with the Dirac point. In fact, special 
measures are required to get them to coincide |21j. 

Whether the condition of weak disorder is realized is 
harder to judge, as very few transport experiments are 
available. It is, however, clear how one would go about 
checking experimentally whether this scenario is actually 
realized. By varying the chemical potential of the surface 
states the maximum of conductance should shift back and 
forth between ^ = and (j) — tt with an amplitude that, 
depending on parameters, is likely to be largest (close 
to e^/h) at the Dirac point (cf. Fig. ^. Alternatively, 
one can change the period of the conductance oscilla- 
tions from (pQ to 00/2 by increasing the surface disorder 
(cf. Fig. [4|. We expect both methods to be practicable 
with current experimental techniques. 

In summary, we have calculated the flux dependence of 
the conductance of disordered topological quantum wires 
in the presence of quantum flux parallel to the wire, mod- 
eling the surface states with a continuous Dirac Hamil- 
tonian. The conductance is found to have the expected 
00 /2 period due to weak antilocalization away from the 
Dirac point and with strong enough disorder such that 
the electron motion is diffusive. A period of 0o is ob- 
tained either at the Dirac point at any disorder strength, 
or away from the Dirac point with weak disorder. While 
in the former case the conductance always has a maxi- 
mum at = 4>q/2, in the latter case the maximum can 



be reached at either </> = or = 0o/2 depending on the 
doping level. We hope that these results will aid further 
transport experiments on the remarkable surface state of 
topological insulators. 

The protection of a mode at 0o/2 is studied using a 
fully 3D simulation in a concurrent work |23| . We ac- 
knowledge conversations with J. Cayssol, D. Carpentier, 
E. Orignac, and A. Vishwanath. This research was sup- 
ported by the U.S. Department of Energy under Con- 
tract No. DE-AC02-05CH11231 (JHB and JEM) and 
the Alexander von Humboldt Foundation (PWB). 
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